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Abstract The pinched/nonpinched classification of intersections of causal singularities of 
propagators in Minkowski space is reconsidered in the context of the theory of asymptotic 
operation as a first step towards extension of the latter to non-Euclidean asymptotic re- 
gimes. A highly visual distribution-theoretic technique of singular wave fronts is taylored 
to the needs of the theory of Feynman diagrams. Besides a simple derivation of the usual 
Landau equations in the case of the conventional singularities, the technique naturally ex- 
tends to other types of singularities e.g. due to linear denominators in non-covariant 
gauges etc. As another application, the results of Euclidean asymptotic operation are ex- 
tended to a class of quasi-Euclidean asymptotic regimes in Minkowski space. 



Introduction 1 

The problem of constructing asymptotic expansions of mul- 
tiloop Feynman diagrams in masses and momenta for various 
asymptotic regimes in a phenomenologically meaningful and 
calculationally useful form is the key analytical problem of ap- 
plied Quantum Field Theory. This is because some sort of as- 
ymptotic expansion seems to always be involved in an impor- 
tant way — explicitly or implicitly — whenever one deals with 
Feynman diagrams in applications (see [1] for a systematic dis- 
cussion). 

The method of asymptotic operation (As-operation for 
short; for a review see [1]) has proved to be the most powerful 
one for solving the asymptotic expansion problem within ap- 
plied QFT. It originated [2]-[5] in the context of studies of 
analytical calculational methods for the coefficient functions of 
short-distance operator product expansion. The resulting algo- 
rithms [3], [6] (extended to arbitrary Euclidean expansions in 
[5]) allow one to obtain contributions to coefficient functions 
from individual diagrams in a maximally simple form, which 
comes about as a natural consequence of a property of perfect 
factorization of the expansions obtained with A.y-operation [3], 
[4]. Due to the perfect factorization, calculations of coefficient 
functions are much simplified and are directly reduced to inte- 
grals for which there exist powerful algebraic algorithms [7]. 
As a result, a class of next-next-to-leading order calculations 
had been reduced to feeding the corresponding diagrams into a 
computer (e.g. [8]-[ll]). A regularization-independent analysis 
[12], [13] demonstrated a remarkable compactness of formal 



proofs within the method of A.s-operation and established the 
latter as a superior alternative to the traditional (BPHZ-type) 
methods in the theory of Feynman diagrams (for a systematic 
dicussion of the differences between the two paradigms see 
[1]). 

It became clear quite early that the method of asymptotic 
operation easily extends to the most general Euclidean asymp- 
totic regimes [4] and that simplicity of the resulting formulas 
remains its characteristic feature [5]. Moreover, the general 
philosophy of the method of A.v-operation [4], [1] is by no 
means specific to Euclidean problems (first examples of non- 
trivial application of the method in essentially non-Euclidean 
situations can be found in [14], [15]). 

The purpose of the present paper is to begin a systematic 
extension of the theory of asymptotic operation to non- 
Euclidean asymptotic regimes. 

Asymptotic operation and causal singularities 1.1 

A necessary preliminary step is to reinterprete the well- 
known pinch/non-pinch classification of intersections of light- 
cone and mass-shell singularities of causal propagators in 
purely distribution-theoretic terms as is necessary within the 
framework of A.s-operation. Indeed, the first question one is 
confronted with is as follows. 

The formulas of As-operation for Euclidean asymptotic re- 
gimes were derived (cf. [4]) using the Wick rotation so that one 
dealt with Euclidean momenta throughout the analytical part of 
the derivation. However, the final formulas are valid irrespec- 
tive of the Wick rotation used at intermediate steps. But then it 
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should be possible to derive the Euclidean Ai-operation di- 
rectly in Minkowski space. 

The technical problem one encounters here is due to the 
more involved structure of singularities of Feynman propaga- 
tors in Minkowski space as compared with the Euclidean 
situation. Indeed, the Euclidean scalar propagator has the form 
(p 2 +m 2 )~ l where p 2 >0 for all p ^ . When formally ex- 
panded in powers of m, this results in singularities at the iso- 
lated point pt 1 = . If there are other integration momenta, 
then the singularity of a single propagator is localized on a lin- 
ear manifold. On the other hand, in a Minkowskian situation 
the standard scalar propagator becomes (- p 2 + m 2 - i0)~ l 
where p 2 can be both positive and negative so that there is a 
singularity localized on the non-linear manifold described by 
the second-order algebraic equation p 2 = m 2 . If m is one of 
the expansion parameters ( m = O(k) ), then the expanded ex- 
pression has singularities on the light cone p 2 = . Moreover, 
the point p = is a singular point of the light cone manifold 
in the sense of differential geometry and should be treated 
separately. The singularity at p = is traditionally referred to 
as soft singularity 1 . Other singular points have the property that 
the propagator near each such point becomes a well-known So- 
hotsky distribution, namely, (x — i0)~ l where x is a properly 
chosen one-dimensional local coordinate. To emphasize their 
analytical nature, such singularities can be referred to as So- 
hotsky-type singularities . To emphasize their physical origin, 
the term causal singularities can be empoyed. While the causal 
- j'O prescription ensures that the propagator itself is a well- 
defined distribution everywhere in the space of p, this is not 
necessarily true for products of such distributions. The causal 
Sohotsky-type singularities may overlap and make the product 
non-integrable by power counting. 

In general, singularities that are non-integrable by power 
counting result in non-trivial contributions to asymptotic ex- 
pansions. Therefore, to obtain a Minkowski-space derivation of 
the Euclidean A.s-operation one has to understand why the in- 
tersections of causal singularities do not contribute in the case 
of Euclidean asymptotic regimes. 

The usual explanation of non-contribution of causal singu- 
larities in Euclidean case is, because the energy integration 
contours are not pinched and one can deform them away from 
such singularities (cf. the derivation of the Landau equations 
that give a criterion of occurence of pinches [16], [17]). But 
this is the same as to say that the Wick rotation can be per- 
formed, and cannot satisfy us. Indeed, eventually one will have 
to answer the question of what happens near the points that 
separate pinched and non-pinched singularities." Then one 
would have to cope with the fact that the latter are disposed of 
via contour deformations while no such trick is allowed for the 
former. To avoid such a situation one has to reinterprete the 
mechanism of non-contribution of non-pinched singularities in 
a manner that avoids contour deformations. 

Moreover, the technique of A.v-operation is, strictly speak- 
ing, not naturally compatible with contour deformations. In- 
deed, whenever distributions are involved one has to deal with 
test functions which in the theory of A.s-operation should be lo- 
calizable — i.e. zero everywhere except a small neighbourhood 



' By association with soft (zero-momentum) photons in QED; analytically, 
this is the hardest singularity generated by a propagator. 
" I first heard this question explicitly stated from G. Sterman. 



of a point which one studies. This, however, cannot be 
achieved with analytic functions, so that analytic continuation 
is no longer allowed. 

Lastly, as was stressed e.g. in [17] (Chapter 2 and refs. 
therein), an accurate derivation of the Landau equations via 
contour deformations in a general case requires a use of the 
notoriously cumbersome methods of algebraic topology. This is 
unsatisfactory because a calculationist would always prefer a 
direct analytical argument at a "microlocal" level — and argu- 
ment that could be directly connected with how one calculates, 
say, asymptotic expansions — to an indirect one such as the 
one based on contour deformations. 

Wave fronts and Feynman diagrams 1.2 

An alternative technique that would allow a direct study of 
overlapping causal singularities is based on the notion of sin- 
gular wave front of a distribution (or simply wave front ) intro- 
duced by Hormander [18] and Sato [19] in the context of stud- 
ies of propagation of singularities of generalized solutions of 
partial differential equations'". Whereas Hormander worked 
with distributions proper [20], Sato considered a related notion 
of hyperfunctions [21]. From a practical point of view, how- 
ever, the two definitions are equivalent (cf. the review [22]). 
The wave front describes singularities of a given function 
(their positions and "orientation", cf. Sec. 4.23 below; for a 
general definition see [22], [23]). The key technical result of 
interest to us here is a criterion of existence of products of dis- 
tributions in terms of their wave fronts (described e.g. in [23], 
Theorem K.45 essentially due Hormander). However, we will 
not need the criterion in full generality, and in Sec. 4 its sim- 
plified version will be taylored to the needs of the theory of 
Feynman diagrams. 

Sato et al. [24] were first to apply the wave-front-related 
techniques to Feynman diagrams. They did so from the point of 
view of the theory of hyperfunctions [21], which was rather 
natural because hyperfunctions can be roughly thought of as 
boundary values of analytic functions." The main result of [24] 
is a description of non-analyticities of Feynman diagrams in 
terms of certain closed formulas (to be compared with the im- 
plicit description via the Landau equations). Despite being in a 
formal sense explicit, such formulas nevertheless seem to be 
less useful in physical applications than, say, the transparent 
interpretation of pinched singularities due to Coleman and 
Norton [26]. v vl ' v " Yet the notion of wave front — if used prop- 
erly, i.e. without excessive generality — is a handy tool in the 
theory of Feynman diagrams due to its visual simplicity. 



"' Cf. the light cone singularities of causal propagators that are special so- 
lutions of the wave equation in coordinate representation. 
1V The relevant physically significant results include e.g. the edge of the 
wedge theorem in the context of studies of dispersion relations (see e.g. 
[25]). 

v I than J. Collins for a discussion of the Coleman and Norton interpreta- 
tion. 

V1 The work of Sato et al. should not be judged too narrowly: they studied a 
new "mechanism of proofs" in an applications-motivated problem in ex- 
actly the same way as physicists would design an experiment to test, say, a 
particle production mechanism — irrespective of its immediate usefulness. 
v " The results of [24] were later reproduced in a series of publications by 
Smirnov (see his compilation [27] and refs. therein) who claimed to have 
used a "different" technique of "singular wave front" of Hormander rather 
than the "analytical wave front" of Sato et al. Such claims cannot be sus- 
tained in view of the exact equivalence of the two versions of the definition 
of wave front in the context of Feynman diagrams, as was emphasized in 
connection with Smirnov' s "results" in [28]. 
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Because the method of A.s-operation iteratively reduces 
studying any however complex geometrical patterns of singu- 
larities to the case of singularities localized at isolated points, 
Sec. 3 contains an exhaustive treatment of one-dimensional So- 
hotsky distributions, the one-dimensional prototypes of causal 
singularities. We consider in detail the issue of when products 
of such singularities exist despite divergence by power count- 
ing. Then this fact is reinterpreted in terms of Fourier trans- 
form thus preparing ground for an extension to many dimen- 
sions. Finally, a simple extension to asymptotic expansions of 
products of Sohotsky distributions is presented. 

Sec. 4 treats the multidimensional case. First, general inter- 
sections of Sohotsky-type singularities are considered, and a 
special version of the highly visual notion of singular wave 
front is deduced in a few easy steps, along with the necessary 
heuristic motivations. A criterion of local existence of such 
products (a special case of the more abstract one due to Hor- 
mander) is derived as an obvious generalization of the one- 
dimensional results, both in an analytical and geometrical 
form. Then it is shown that the same criterion holds for the 
case of osculating singularities; at this point the heuristic 
power of the notion of wave front becomes apparent. 

Finally, in Sec. 5 we turn to multiloop diagrams and present 
a one-step derivation of the Landau equations. Then we briefly 
discuss how the developed formalism helps to understand the 
structure of singularities in the theory of non-covariant gauges, 
and how the results of Euclidean theory of asymptotic opera- 
tion are extended to the class of quasi-Euclidean asymptotic 
regimes. 

Causal singularities and asymptotic 
operation 2 

Notations for integrands 2.1 



Purpose 1.3 

It was already emphasized elsewhere [1] that the mathe- 
matical apparatus of asymptotic operation includes a collection 
of techniques that cannot be found in any single source on dis- 
tribution theory. 1 This is why it would be useful to present an 
account of the corresponding techniques in a form taylored to 
the specific needs of theorists that have to deal with Feynman 
diagrams in applied problems. The goal of the present paper 
from a technical point of view is to do that with respect to the 
technique of wave fronts. 

In this respect note, first, that applications to Feynman dia- 
grams do not require the use of wave fronts in full generality. 
Second, whereas the usual treatments consider "static" prod- 
ucts of singular functions, we wish to study asymptotic expan- 
sions so that not only existence of products has to be estab- 
lished but also the fact that they possess the required asymp- 
totic approximation properties. All the necessary mathematical 
apparatus will be explained in detail. 

From the point of view of the theory of multiloop diagrams, 
we are going to: rederive the Landau equations [16], [17] in a 
purely distribution-theoretic fashion, i.e. avoiding the use of 
contour deformations as is appropriate within the theory of As- 
operation; demonstrate how the technique can be applied to 
analysis of singular functions of a more general form than the 
standard causal propagators (cf. the problem of singularities in 
non-covariant gauges [31]); extend the Euclidean theory of as- 
ymptotic operation [4], [5] to the simplest class of non- 
Euclidean asymptotic regimes (the so-called quasi-Euclidean 
ones). 

It should be emphasized that the motivations of this work go 
somewhat beyond providing a "better derivation" of the Lan- 
dau equations. The problem of asymptotic expansions of mul- 
tiloop diagrams in non-Euclidean regimes is a notoriously 
complex one, and complex problems are solved by building an 
adequate intuition based on — and together with — and ade- 
quate formalism. So our aim is to present a more flexible tech- 
nique than the one based on contour deformations and to re- 
consider the classification of singularities in a manner suitable 
for a subsequent construction of non-Euclidean infrared R- and 
A.v-operations. The fact that the Landau equations reemerge as 
a result is in principle of little consequence. 

Plan 1.4 

The plan of the paper is as follows. In the descriptive Sec. 2 
we establish the context in which to study non-Euclidean As- 
operation. First the systematic notations are summarized to de- 
scribe the products of singular functions we wish to deal with. 
Then a summary of A.s-operation is presented emphasizing the 
point that its structure is independent of whether one works in 
Euclidean or Minkowski space. Further, we identify the non- 
linear causal singularities and briefly discuss the distinctively 
non-Euclidean geometrical phenomenon of osculating singu- 
larities. The conventional treatments avoid discussing osculat- 
ing singularities at a microlocal level altogether, but the issue 
cannot be circumvented if one aims (as we do in the theory of 
A.s-operation) at a complete solution of the asymptotic expan- 
sion problem. 



' The original source on the theory of distributions is [20]. Good textbooks 
are [29] (with many examples to different physical problems) and [30] 
(with proofs of the general results of the theory of distributions that avoid 
the abstract theory of topological spaces). Wave fronts are briefly dis- 
cussed in [23]. 



Our notations will generally agree with [4]. Let us briefly 
summarize them together with whatever extensions are needed 
in the non-Euclidean case. This will also introduce one into the 
context of the problem. 

The object one deals with is G(p) , the integrand of a mul- 
tiloop diagram G, whereas p = (p l ,...p l ) is the collection of 

all its integration (loop) momenta: dp = A D p x d D p 2 ... . The 

space of p is denoted as P. The concrete value of D (the num- 
ber of scalar components of each p i ) is not important. It is of- 
ten convenient to ignore the structure of P altogether and to 
consider it as an abstract vector space of finite dimensionality. 
The diagrammatic interpretation of G will not be important 
either and we will treat G(p) simply as a product of singular 
factors. Its structure is described below. 

A general scheme is that one starts with G(p) and con- 
structs a distribution from it. The conditions governing such a 
construction are always local. This means one should first per- 
form such a construction for all sufficiently small O 
(neighbourhoods of singular points). After that the transition to 
a distribution defined on the entire P is achieved via the stan- 
dard trick of the decomposition of unit [29], [30], For this rea- 
son, we focus attention only on what happens in a small open 
region O. 

Furthermore, we will ignore the polynomial factors in the 
numerator of G(p) that are due to interaction vertices or non- 
scalar particles. This is because the A.s-operation, by definition, 
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constructs a complete asymptotic expansion in powers and 
logarithms of the expansion parameter. Such expansions are 
unique and, therefore, commute with multiplications by poly- 
nomials [4]. This is a general property that remains valid in the 
non-Euclidean case. Therefore, the presence of polynomial 
factors is of no significance as far as the analytical mechanism 
of As-operation is concerned. They play no role in the study of 
the geometry of causal singularities and will be ignored. 

So, the products of singular factors that we consider in this 
paper have the following form 

G 0>) = n seG sO>). 2.2 
where 

g(p)=A g (p) = (D g (p)-iO)~" g . 2.3 

D g (p) is whatever expression happens to occur in the de- 
nominator of the g-th factor. Eq.2.3 fixes the relative sign of 
D g (p) with respect to the infinitesimal imaginary part. In this 

work we consider only such regions O C P that the change of 
sign does not occur within O in all factors in the product. 

Structure of denominators 2.4 

In the case of the standard scalar propagator, 
D g (p) = -l 2 g (p)+m 2 g , 2.5 

where m g may or may not be zero. The quantity l g (p) is the 

momentum flowing through the g-th line of G; it is a linear 
combination of p i and the external momenta. The concrete 

form of I is determined by the topology of the diagram and 

the pattern of its external momenta. When studying the general 
analytical mechanisms, we make assumptions neither about the 
topology, nor about the pattern of external momenta (i.e. 
whether they are put on mass shell, taken at zero values etc.). 

The construction of A.s-operation requires to perform vari- 
ous expansions of the factors [1], [4]. This explains why we 
allow the powers n g to differ from 1. Moreover, in the non- 
Euclidean case, such expansions may change the form of de- 
nominators. In particular, there can emerge denominators with 
a linear dependence on the momentum, 

D g (p) = n g p + m 2 g , 2.6 

where n g is a parameter vector, etc. The examples where such 

factors occur are: the string operators in the context of the Su- 
dakov problem etc. [32]; the propagator of the heavy quark in 
the effective heavy quark theory (cf. e.g. the effective QCD La- 
grangian in [33]). The non-covariant gauges [34] also involve 
factors with denominators similar to 2.6 but with the sign of 
the imaginary part different in different regions of P. 

hi view of the variety of possible denominators, it makes 
sense to keep irrelevant details out of the way and simply re- 
gard D g (p) as a smooth function of p. 

The expansion parameter is denoted as k. The dependence 
on/c of various quantities is indicated, as usual [4], with an ad- 
ditional argument, e.g. G(p,/c) , D g (p,K) etc. Its particular 

form is determined by the specific problem and the asymptotic 
regime chosen. 



General structure of the asymptotic operation 2.7 



It is important to understand that the general structure of the 
A.s-operation is the same irrespective of whether one deals with 
Euclidean or non-Euclidean problems. This is because neither 
the motivations for the distribution-theoretic point of view pre- 
sented in Sec. 4 of [4], nor the extension principle of Sec. 6 of 
[4], nor even the general structure of the arguments of Sec. 11 
of [4] which established the structure of A.v-operation, depend 
on that assumption. Let us summarize the prescriptions of As- 
operation emphasizing this fact. 

The asymptotic operation As K (with respect to a concrete 

small parameter k which implies a concrete choice of the as- 
ymptotic regime) is meant to be applied to a product of singu- 
lar factors G(p,K) regarded as a distribution in the integration 
momenta p. It is then supposed to yield an asymptotic expan- 
sion of G(p,K) in powers and logarithms 1 of k in the sense of 
distributions. At an operational level, the construction involves 
the following steps: 

(i) Applying the formal Taylor expansion T K in powers of k to 
each factor A g (p,K) in the product, 

Govc)=n seG s(p^) 

- t * ° n geG *o»."> = n geG t , ° s(p^) , 2.8 

which results in a formal series in powers of k, identification of 
potentially non-integrable singularities in the resulting expres- 
sion and their geometrical classification. Here one determines 
the number and geometric structure of the manifolds in P on 
which those singularities are localized. One here deals with an 
hierarchy of intersecting manifolds of various dimensionalities. 
Note that all such manifolds occur as intersections of singular 
manifolds for individual factors g(p,K) . 

(ii) Determining which factors from the product contribute at 
each singular point. Here one establishes a correspondence 
between the singular manifolds n y , the subproducts y C G 

consisting of factors that are singular there, and the diagram- 
matic interpretation of such subproducts as some kind of sub- 
graphs of the initial Feynman graph." 

(iii) Studying the analytical structure of the singularities (some 
form of power counting etc.) and constructing a subtraction 
procedure that transforms the singular and, in general, non- 
integrable formal expansion into a correctly defined distribu- 
tion on the entire P. 1 " As in [13], we denote it as R : 



1 The powers need not be integer, e.g. square roots are allowed etc. We are 
talking about "powers and logarithms" although the analytical fact that the 
expansions we deal with run in powers and logarithms of the expansion pa- 
rameter and do not contain, say, log log K terms, cannot be decided a 

priori. However, such an assumption is not used in the actual construction 
of As'-operation and its correctness is established along the way. Note also 
that proving the power-and-log structure of expansions is a much simpler 
problem than establishing their precise diagrammatic form. References to 
some such studies can be found e.g. in [4]. 

" Note that the resulting diagrammatic images are determined by the un- 
derlying analytical structures (in particular, the asymptotic regime) and 
need not fit exactly into the standard categories of the theory of graphs. In 
fact, a graph-theoretic characterization, even if possible, may be too cum- 
bersome to be useful, or simply irrelevant (except in a very special context, 
e.g. when rewriting known results in a "rigorous" form; cf. [27]). 
"' In certain cases one may do without such a subtraction procedure if a 
convenient regularization is available (the dimensional regularization is 
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n seG To ^)-R°n, eG To ^) 29 

The subtraction procedure R should be minimal in that it 
should not modify the leading power behaviour near singular 
points. It is analytically somewhat similar to the Bogoliubov R- 
operation in coordinate representation. 

(iv) Constructing the counterterms E y (p,/c) localized at the 
singular manifolds n y (thus fixing the finite arbitrariness in 

the definition of R ): 



^ ° ,efi ^ ° 8(P> K ) 



-*R° l\T K og( p , K )xE Y (p, K ) = As K oYlg( p , K ). 2.10 

g EG\y ' gEG 

The inclusion of E y (p,/c) is necessary in order to transform 

the expression into a correct expansion in the sense of distri- 
butions. E y (p,/c) contain all the non-analytic (logarithmic, 

square root etc.) dependence on k. They are localized on the 
singular manifolds 7t and, therefore, are linear combinations 

of (5-functions. The number of derivatives of the (5-functions is 
determined by the power counting at the previous step, 
whereas the non-trivial coefficients are found using the con- 
sistency conditions [4]. The coefficients of the counterterms of 
the A.v-operation are detemined uniquely [4]. 

Note that one needs some form of explicit parametrization 
of n y to write down the structure of (5-functions (their argu- 
ments etc.). 



Singular manifolds 



2.11 



Strictly speaking, in the context of asymptotic expansions 
one deals with the singularities of the formal expansions of the 
product G(p,/c) in k. Note, however, that the structure of 
those singularities is the same (up to the strength of singulari- 
ties) for all terms of the formal expansion. This is because the 
singularities are due to denominators, and the terms resulting 
from expansion of a factor have the same expression raised to 
different powers in the denominator. Therefore, it is conven- 
ient first to ignore the "dynamic" aspect of the problem (i.e. 
the fact that one deals with expansions) and consider the 
problem as "static" (i.e. existence of products only). After that 
we will see what modifications are needed to take into account 
approximation properties etc. 

The singularities of the product G(p) are due to the fact 
that each of the denominators may have zeros, and the singula- 
rities of individual factors in the space P may overlap in a non- 
trivial fashion. One wishes to study such overlapping singula- 
rities based on the information about the individual factors. 



the preferred choice). Then the task of recasting the final answers into an 
explicitly convergent form is postponed till the end of the construction — 
cf. [4], [5] where such an approach was motivated by calculational con- 
venience of final formulas in a situation where the subtraction procedure is 
understood well enough from analogy with the UV renormalization in co- 
ordinate representation (cf. [36], [12]). Such an approach, however, may 
not be always desirable, especially in the non-Euclidean problems where 
the multi-level recursion structure of the problem with secondary expan- 
sions obscures the connection between various types of divergences that 
occur at intermediate steps. Moreover, the dimensional regularization may 
not always work in non-Euclidean situations [14]. 



The singularities of the g-th factor are localized on the mani- 
fold n g C P described by the equation D g (p) = 0: 



it ={p\DJp) = 0} . 



2.12 



The singularity corresponding to the g-th factor is referred to 
as elementary singularity . 

Note that Ji g may be a smooth manifold everywhere in P 

(e.g. the cases 2.5 with m g ^ and 2.6). On the other hand, 

the points where l g (p) = for 2.5 with m g = (the so-called 

soft singularity) correspond to the apex of the light cone where 
the light cone is not a smooth manifold. More generally, such 
special singularities occur whenever the gradient nullifies, i.e. 
VD (p) = , where the gradient operator is 



2.13 



_a a_ 

dpi dP2 

If D g (p Q ) = but VD g (p Q ) ^ , then there is a small neigh- 
bourhood O of p Q in which the local coordinates 
p = p(x,y,...) may be chosen so that D g (p(x, y, ...)) = x . 
Then the factor becomes 



A g (p(x,y,...)) = (x-i0) n s . 



2.14 



We call the distributions that can be locally reduced to such 
form Sohotsky-type distributions . Their integrals with test 
functions are well-defined despite the fact that they are diver- 
gent by power counting (see also Sec. 3). 

hi this paper we consider only intersections of causal sin- 
gularities. Formally, this means that all Ji g are smooth mani- 
folds within O. A sufficient condition to ensure that is to re- 
quire that 

VD g (p)*0, Vp£(9, Vg . 2.15 

This means that we exclude from consideration the soft singu- 
larities — their study requires additional techniques that would 
lead us outside the scope of the present paper. The restriction 
will remain valid throughout this paper. 



Conventions for the gradient V = d/dp 



2.16 



Geometric representation of the gradient vector V£)(p) in 
the space P of independent variable p depends on how one 
chooses the scalar product in P. This is because the gradient is, 
strictly speaking, a vector from the co-space P' whose embed- 
ding into P is a matter of convention. However, the properties 
we are going to consider are expressed via relations that are 
linear in V (e.g. VD(p) ^0) and are, therefore, unaffected by 
the conventions adopted. This allows some flexibility of nota- 
tion depending on specific purposes: To build insight into the 
nature of intersections of Sohotski-type distributions, we will 
consider simple examples, in which it is convenient to regard 
P is a Euclidean space. Then the gradient vector VD(prj) l& 
normal to the singular manifold D(p) = 0. But after the gen- 
eral theorems are understood and we turn to concrete examples 
in the context of Feynman diagrams and Minkowski space, it is 
convenient to proceed in a formal algebraic manner, hi par- 
ticular, p is then an array of Minkowskian vectors (the collec- 
tion of loop momenta) and is of a mixed Euclidean/Minkowski 
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type, e.g. p- p' ' = Y,i Pi Pi fi ■ Then the gradient VD(p ) may 

happen to be tangential to the corresponding manifold 
D{p) = (for instance, if p is a Minkowskian 4-vector and 

D(p) = p 2 then VD(/? ) = 2p(^). But that is a purely formal 
fact of no analytical consequence. 

Singular manifolds and singular subgraphs 2.17 

Consider an arbitrary subset of lines, yCG, and define: 



Example: self-energy 



2.21 



gEyCG 



In other words, n y is the collection of points through which 

pass the singularities of each factor g £ y . It is called singular 
manifold corresponding to the subproduct y . 

As in the Euclidean case, one may have ji y = n , even 

when y ^y' . However, the analytical nature of the singularity 
of the entire product G(p) on n y is determined by all the 

elementary singularities ji g Dji y and by all factors g(p) that 

are singular on the entire n y , i.e. ji g D n y . So, whenever one 

deals with a singular manifold Ji y , the corresponding y C G 

contains all g such that Ji g D Ji y . The subproduct 

y(p)=f] gCy g(p) 2.19 

that contains all the factors from G(p) that are singular on 
ji y is described as complete . The singular manifold Ji y to- 
gether with the associated complete subproduct y will be called 
singular subgraph although at this stage we do not associate 
any graphical image with it. 

The collection of all singular subgraphs of G is denoted as 
S[G]. 



Osculating singularities 



2.20 



In the Euclidean case, there were two basic geometric pat- 
terns of intersections of elementary singularities: factorizable 
and non-factorizable. Because now singular manifolds are non- 
linear, an additional pattern has to be dealt with, namely, os- 
culating 1 singularities. 

Let g(p) be an intersection point of several non-linear 

manifolds ji g , g £ y , and let O be its small neighbourhood. 

One can always find a smooth deformation of coordinates 
p -* p' within O that transforms any one of the singular mani- 
folds Ji g into a linear manifold so that the corresponding fac- 
tor in the product takes the form (y g p' — iO) " g where v g is a 
constant vector. There are two cases possible: 

• A general case is when all the factors g(p)E.y that are 
singular at p Q can be simultaneously linearized. 

• When a simultaneous linearization is not possible, one deals 
with osculating singularities . 



Consider the integrand of a one-loop scalar self-energy in- 
sertion: 

i i 

2.22 



-p z +mf-i0 -(p-qf +m$-i0 

where m i ^ , p is a 4-dimensional integration momentum 
while q is the external momentum. The integral of 2.22 over p 
(after appropriate UV renormalization) is a function of q 2 and 



213 the masses. The interesting asymptotic regimes are as follows: 

(i) q 2 -» - oo (the Euclidean regime), 

(ii) q 2 -» +oo (the quasi-Euclidean regime; cf. Sec. 5.9); 

(iii) q 2 ->(m l +m 2 ) 2 (the threshold regime). 

To study (i), one chooses q to be purely space-like and then 
Wick-rotates p into the Euclidean region. Then the formalism 
of [4] becomes applicable. 

The case (ii) is more interesting because from the point of 
view of momentum space, the picture here is much different 
from the Euclidean case. Indeed, choose q to be purely time- 
like, q = (^q 2 ,0,0,0 j . An immediate Wick rotation is prohib- 
ited. Rescale the dimensional parameters, p -» p(q 2 ) 112 and 
mf ->mfq 2 , so that the regime (ii) is replaced by m i -* for 

q fixed, with both masses going to zero at the same asymptotic 
rate. Within the framework of the theory of As-operation, one 
constructs the corresponding expansion of 2.22 in the sense of 
distributions. The first step of the construction is to identify the 
singularities of the formal Taylor expansion, whose terms have 
the following general form: 

1 1 



(- P 2 -;or> (-(p-qV-iOy 



2.23 



The singularities are localized on two light cones, one shifted 
with respect to the other by q as shown in 2.24: 



p = q 




P = 



2.24 



1 i.e. "kissing" (from Latin). 



Two counterterms will be needed for the singularities at the 
apices of the light cones. Recall that in the Euclidean regime 
(i) one would have to deal with these two points only. In the 
present case, however, there are also singularities localized on 
the light cones away from their apices. Thus, in a small neigh- 
bourhood O of a general point A of any of the light cones the 
local coordinates can be chosen so that the corresponding factor 

becomes (x — i0)~ N ' , which is non-integrable by power 
counting for N t large enough. The prescriptions of As- 

operation would require counterterms to be introduced for such 
singularities. However, it turns out that such light-cone singu- 
larities do not contribute to the final answer. The usual argu- 
ments explaining this rely on deformations of the contour of 
integration over p° into the complex plane. But for reasons 



6 



F. V. Tkachov: Landau equations.. 



discussed in the Introduction we would like to find an expla- 
nation for the absence of the counterterms in such cases di- 
rectly in real Minkowskian momentum space. 

The threshold regime (iii) is the most difficult one. Indeed, 
the small expansion parameter k can be introduced by repre- 
senting q as follows: 



q = (mj +7W2 + /c, 0,0,0) 



2.25 



Then the expansion in k -» is equivalent to the expansion in 
q 2 -> (m l +m 2 ) 2 . 

A typical term of the formal expansion of 2.22 in k has the 
following structure of denominators that determine the struc- 
ture of singularities: 

1 



[-2m 1 r -(r ) 2 +/- 2 -i0] 
1 



l+lrr^r" - (r u ) z + r z - i0]" 
where we have introduced the variable r: 



P = (P°,P)-- 



(m, +r°,r) 



2.26 



2.27 



The singularities of the two factors are localized on two hyper- 
boloids which touch (osculate) at the point r = 

(p = q=ql =n )- 



p 2 = mf, p°>0 




2.28 

By analogy with the preceding example one has to conclude 
that the general points of the hyperboloids need not be supplied 
with counterterms when the expansion in the sense of distribu- 
tions is constructed. But the argument does not work for the 
point r = (p = q) . On the other hand, it is well-known from 

studies of analyticity properties of Feynman diagrams that the 
asymptotic regime under consideration corresponds to the point 
of non-analyticity; therefore, the expansion must contain a non- 
analytic term. In view of the fact that non-analytic dependences 
are contained in the coefficients of counterterms of the As- 
operation (cf. [4]), one concludes that it is exactly the coun- 
terterm for the singularity at r = (p = q) that should contain 

the non-analytic threshold behaviour. 

Consider, however, the expansion in k defined as follows: 



-«2+/c, 0,0,0) . 



2.29 



Then there will also be a point where the two hyperboloids 
touch (cf. 2.30; only the osculating upper hyperboloids are 
shown; the dotted line shows the position of the unshifted hy- 
perboloid ). 



But according to the known analyticity properties, the expan- 
sion in k near 2.25 should not contain non-anaytic terms. This 
means that, unlike the preceding case, the osculation point of 
the two hyperboloids does not require a non-trivial coun- 
terterm. 



Example: collinear singularities 



2.31 



Another rather typical geometrical pattern of intersection of 
causal singularities occurs in the following product of two 
massless propagators: 



1 



1 



[-p 2 -i0fi [-(p-q) 2 -i0fi 



2.32 



where q = (both p and q are 4-dimensional Minkowskian 
vectors). Such a product may have resulted from a formal ex- 
pansion in one of the high-energy regimes (cf. [35]). The sin- 
gularities of the two factors are localized on the two light cones 
as shown in 2.33. 




2.33 



The first factor is singular on the light cone p 2 = while the 
second, on a similar light cone whose apex is shifted to the 
point p = q . The two light cones intersect over the line 



L = {p\p = zq] 



2.34 



(i.e. p should be collinear to the light-like momentum q which 
is the reason why this is called collinear singularity ). Note 
that, in general, an intersection of two 3-dimensional manifolds 
is 2-dimensional, while in the present example it is 1- 
dimensional. To see what happens here, consider a cross sec- 
tion of the space of q with a hyperplane p° = const > q° 
(represented by the line AA in 2.33) as in 2.35. 



= 




(P-q) =0 



2.35 



The intersections of the two light cones with this hyperplane 
are spheres that osculate at a point. When the cross section 
shifts along p° , the spheres inflate/deflate and shift, but there 
is always one osculation point. If < p° < q° then the two 
spheres are side-by-side, otherwise one is inside the other; the 
two cases correspond to pinched and non-pinched singularity. 
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One-dimensional case 3 

Within the formalism of asymptotic operation, all problems 
are recursively reduced to investigation of one-dimensional 
singularities. Therefore, it is convenient to review the proper- 
ties of the relevant one-dimensional distributions and remind 
ourselves of some basic properties of Fourier transforms that 
will be needed in the discussion of the multidimensional case. 
The Sohotsky distribution 3.1 

Let p be a one-dimensional integration variable running 
from —oo and +00 . Recall the familiar one-dimensional So- 
hotsky distribution 

' lim 



p — iO tj^O p — if] 



3.2 



which is well-defined despite the fact that the singularity is 
logarithmically divergent by power counting. As usual, the 
limit should be understood in the sense of the distribution the- 
ory, i.e. the above expression is not expected to be interpreted 
otherwise than in integrals with test functions <p(p) (that are 

smooth everywhere and non-zero only in a compact region): 

+ 00 +00 

1 1 
f dp <p(p) * lim [dp <p(p) . 3.3 

J p — lV n^O J V — m 

— 00 —00 

Note that both the integral and the limit on the r.h.s. exist for 
any such <p(p) . 

The well-known Sohotsky formula reads: 

— ! — =PV— + ijzd(p) , 3.4 
p — iO p 

where the first term is defined with the usual principal value 
prescription. 



Products of Sohotsky distributions 



3.5 



The next point is existence and non-existence of products of 
Sohotsky distributions. For instance, the product of the Sohot- 
sky distribution and its complex conjugate, 



1 1 

-X- 



3.6 



p — iO p + (0 

does not exist as a distribution defined on all test functions 
unlike the square 



1 1 

-X- 



p — iO p — i'O 



3.7 



which does. 

In the context of analytical properties of Feynman diagrams 
(cf. e.g. [17]), this phenomenon is well-known and described in 
terms of pinched (Eq. 3.6) vs non-pinched (Eq.3.7) integration 
contours. Indeed, if one integrates Eq.3.7 with a function that 
is analytical in a neighbourhood of p = , then the Cauchy 
theorem allows one to simply deform the contour of integration 
into the upper complex halfplane near p = because the poles 
of the two factors are on the same side of the contour; then the 
integrand will have no singularities at all on the integration 
contour. In the case of 3.6 this is not possible — the contour is 
"pinched" by the two poles. 



Interpretation in terms of Fourier transforms 



3.8 



Following, essentially, Hormander (see [23] and refs. there- 
in), consider the products 3.6 and 3.7 from the point of view of 
Fourier transformation. Then existence of products of functions 



that are singular at p = is determined by the behaviour of 
Fourier transforms at infinity. Such a behaviour may be differ- 
ent in different directions, and under certain conditions this 
may result in a well-defined product even if the factors are sin- 
gular. 

First recall that Fourier transforms of test functions are 
smooth functions that exhibit rapid decrease in any direction 
(i.e. decrease faster than any inverse power of their argument). 
However, Fourier transforms of distributions need not be de- 
creasing. On the other hand, they may behave badly in some 
directions, and exhibit a rapid decrease in others. For instance, 
by straightforward integration one sees that 



1 + °° 

— = i jdxe- ikx e-° x e(x>0) 



p-iO 



3.9 



• +00 



The Fourier transform 9(x > 0) remains constant as x - 
but decreases rapidly as x -» - 00 . 

For the complex-conjugate distribution the pattern of dan- 
gerous and non-dangerous directions is mirrored. 

Note that raising the l.h.s. of 3.9 to an arbitrary power 
n > results in a power of x on the r.h.s.: 



1 



p-iO 



\n +00 

= ( jdx 



e i P x (-ix) n ~ l e- 0x 6(x>0) 



3.10 



The important point is that the growth of the integrand corres- 
ponding to a power singularity is no more than polynomial. 

Further, consider the integral of 3.7 against a test function 
(p{p) . Using 3.9, we find: 



1 1 

[dp -X -<p(p) 

J p-iO p-iO 



p — iv p- 

+00 +00 



= const jdx jdx' 0(x>0) Q(x >0) y(x + x) . 



3.11 



-00 —00 



Thus, the convergence of the integral on the l.h.s. around 
p = has transformed into convergence of the integral on the 
r.h.s. at infinity in various directions. (Since the Fourier trans- 
form (p(x) is a smooth function, integration over finite regions 
is not dangerous.) One sees that the ^-functions ensure conver- 
gence in all directions other than those in the positive quadrant 
x > 0, x' > . However, convergence in those directions is en- 
sured by the properties of (p(x) , as explained above. 

One concludes that the integral 3.11 exists for any test 
function, therefore, the product 3.7 is a well-defined distribu- 
tion. 

On the other hand, a similar calculation for the product 3.6 
gives: 

+ 00 



fdp— ^-X— J— p(p) 

J p-iO p + iO 

-00 

+ 00 +00 

c jdx jdx' 6(x>Q)6(x' <Q) &(x + x') . 



3.12 



-00 — 00 



One can see that there are directions in the plane (x,x : ) — 
e.g. x = A + x Q , x' = — A + Xq — in which the ^-functions 
are both equal to 1 while ~p{—x + x') = (f(—x Q + Xq) = const as 
A -» +00 . The integral on the r.h.s., therefore, cannot be con- 
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vergent for all test functions, and the product 3.6 is not a well- 
defined distribution. 

As a further excercise, replace the second factor by, say, 
d(p) . Then the second 0-function on the r.h.s. of 3.11 is re- 
placed by a constant, and again, the integral does not exist for 
some test functions: the product of a (5-function and a Sohot- 
sky-type distribution is ill-defined. 

On the other hand, if one of the distributions is a function 
that has any number of derivatives in a small neighbourhood 
around p = , then there are no dangerous directions associ- 
ated with it at that point, so that its product with any distribu- 
tion exists around p = . 

"Natural" existence 3.13 
It is important to realize that any definition of an integral in 
infinite bounds contains — explicitly or implicitly — some sort 
of a cutoff, say, \x\< A, A -* oo . Such a cutoff plays the role of 

a natural regularization. Then the absolute convergence of the 
integral 3.11 (established e.g. by power counting) implies that 
whatever (correct) regularization method one uses to give 
meaning to the product of the two Sohotsky distributions on the 
l.h.s. of 3.11, the limiting value for the integral on the l.h.s. 
will exist for any test function and it will be the same for all 
regularizations. This gives a precise meaning to the phrase that 
the product in such a case is "naturally defined" in the sense of 
distributions. 



Asymptotic expansions 



3.14 



A single distribution 

As a warm-up excercise, consider the distribution 
1 



p+m — iO 



3.15 



3.16 



where m is a real parameter. Its expansion in the sense of dis- 
tributions as m -» is: 



Eq.3.16 - 2 



{-m) n 



B+l 



3.17 



Proof. It is sufficient to consider only the case n = . As 
usual, the expansion is sought in the form 



1 



_ 1 



p+m — iO p — iO 



+ c(m)d(p)+o(\) 



3.18 



and the coefficient c(m) is determined from a consistency 
condition (in a manner similar to the Euclidean case; see 
sec. 7.4 of [4]): 



c(m) = j dp 



1 



1 



p+m — iO p — iO 



+ o(l) = + o(l) . 3.19 



The integration can be performed using the explicit represen- 
tation 3.4 and symmetric cutoffs, \p\< A . 

The expansion 3.17 can be regarded as the Taylor expan- 
sion of a translation along the p-axis. This explains the ab- 
sence of non-trivial counterterms with logarithmic dependence 
on m : an infinitesimal translation involves differentiation, and 
differentiations in p of the distribution (p±iO)~ n are per- 
formed in a straightforward manner. 

Two factors 3.20 
Consider a product of two factors like 3.16 with their sin- 
gularities separated by a distance o(m) . Then one has the fol- 
lowing expansion: 



- 1 



(-m) n 

-.n+2 



p-iO p + m-iO m ^ ^ (p-iO) 



3.21 



It can be proved in a way similar to the preceding example by 
using the procedure described in sec. 7 of [4]. Another way is 
to use partial fractioning, 



1 



1 



_ 1 



1 



1 



3.22 



p — iO p+m — iO myp — iO p+m — iOj' 

and reduce this example to the case of one factor already con- 
sidered. 

Note that if one changes the sign of the imaginary parts of 
one of the factors, then non-trivial m -dependent counterterms 
appear in the expansion. For example: 

(-m)" 

m 



1 



1 



-2 



p + iO p + m-iO m ^ ^ (p-iO) 



, n +2 



3.23 



The simplest way to obtain this result is via a trick — switch- 
ing the sign of imaginary part of the first factor using the So- 
hotsky formula 3.4. However, one can also derive it in a sys- 
tematic way as in the preceding examples. 

The conclusion is that the expansion 3.21 looks as if obtain- 
ed by simply taking a product of the first factor (whose expan- 
sion in m is trivial) and the expansion of the second factor, 
Eq.3.17. Each term of such a product is a well-defined distri- 
bution. However, strictly speaking, existence of such a product 
does not guarantee that it will possess the necessary approxi- 
mation properties to be the asymptotic expansion of the initial 
m -dependent product (the l.h.s. of 3.21). But it turns out that it 
does. Let us understand why this is so using the techniques of 
Sec. 3.8: in a general multidimensional situation, explicit cal- 
culations similar to those above may not be possible. 



Interpretation in terms of Fourier transforms 
Similarly to 3.11, one has: 

+ oo 

jdp^-^X-^ 

— 00 

+ 00 +00 

oc jdx j dx'8(x>0)[8(x'>0)e~ imx ']^(x + x') . 



3.24 



3.25 



-00 — 00 



If one performs the Taylor expansion in m on the r.h.s., the 
factor exp(-i'mx') (which is bounded by 1 everywhere) trans- 
forms into a polynomial of x' . The remainder term is then 
bounded by something like 0(m N x' N ) which exhibits a poly- 
nomial growth as x' -» ±°° — but not faster than polynomial. 
Therefore, the pattern of behaviour at infinity (polynomial 
growth vs rapid decrease) is the same both for non-expanded 
factors, the expanded factors, and the remainder term. This al- 
lows termwise integration of the product of expansions with 
preservation of approximation properties. 

One sees that the "natural" existence (in the sense of 
Sec. 3. 13) of the product of expansions of individual factors 
implies that such a product is indeed a correct asymptotic ex- 
pansion for the initial non-expanded product; in other words, 
taking product and performing expansion are two commuting 
operations in such a case. 

The above examples provide enough motivation for analog- 
ous multidimensional constructions to which we now turn. 
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Products of Sohotsky distributions 

in d> 1 



1 +0 ° 
= ; fdlf 

L — iO J 



-ilA-OA 



4.8 



We are going to generalize the notion of dangerous direc- 
tions to the multidimensional case. We first consider the case 
of linear denominators (Sec. 4.1) and derive a simple geometric 
criterion of existence of products of Sohotsky-type dis- 
tributions. We then consider smooth deformations of coordi- 
nates (Sec. 4. 14) and see what form the criterion takes. The 
notion of singular wave front emerges naturally (Sec. 4.23) and 
allows one to easily obtain a criterion (essentially due to Hor- 
mander; cf. Theorem IX.45 in [23]) for the general non-linear 
case including the so-called osculating singularities (Sec. 4.34). 
Finally, the extension of the criterion to the case of asymptotic 
expansions of Sohotsky-type distributions becomes straight- 
forward (Sec. 4.41). 



Linear case 



4.1 



Two factors 



4.2 



Consider a product of two Sohotsky-type distributions inte- 
grated against a test function <p( p) localized in the region O: 



J" dp <p(p) 



1 



(vjp - (0)(v 2 p - (0) 



4.3 



For simplicity of notation, we do not consider general powers 
explicitly. The distributions with -HO can be taken into ac- 
count by changing signs of (some of) v,- . 

One can see that the product is a well-defined distribution 
unless the two vectors v x and v 2 look in opposite directions. 

Indeed, if they are linearly dependent then the problem effec- 
tively degenerates into a 1-dimensional one: if they look in the 
same direction then the product is a well-defined distribution; 
if they look in opposite directions then the product is non- 
integrable. If the two vectors are linearly independent then one 
deals with a direct product of two distributions that are well- 
defined separately, and such a product is a well-defined distri- 
bution too. 



Three factors 



4.4 



Under the same assumptions, consider the product of three 
Sohotsky-type distributions: 



J" dp <p(p) 



1 



(vjp - i0)(v 2 p - i0)(v 3 p - iO) 



4.5 



The singularities of the g-th factor are localized on the linear 
manifold 



4.6 



First, without loss of generality one assumes that the origin is 
the only point of intersection of the three manifolds. Second, 
one has to assume that the product of each pair of the three 
factors is well-defined as a distribution everywhere. A standard 
reasoning [4] shows that the product 4.5 is then integrable with 
any test function that is equal to zero in any small neighbour- 
hood of the origin of coordinates. It remains to study the be- 
haviour of the product near that isolated point. 



Existence of the product: analytical criterion 

Using the representation 



4.7 



one can rewrite 4.5 as follows: 

+ 00 +oo +oo 

Eq.4.5 =i 3 JdAj jdX 2 JdA 3 £(v x Aj +v 2 A 2 + v 3 X 3 ) , 4.9 
ooo 

where !p(x) is the Fourier transform of <p(p) : 

4.10 



£(x)ocJdpe-*"V(/>) 



One can see that the rapid decrease of ~p makes the integral 
convergent at large X g in any directions if the following con- 
dition is satisfied 

v x Aj + v 2 X 2 +V3A3 * 0, forallA g >0. 4.11 

Then the product in the integrand of 4.5 exists in the sense of 
distributions. 

Note that for, say, X 3 = , the first inequality in 4.11 de- 
generates into the condition of existence of the product of the 
first two factors. Therefore, instead of making the assumption 
of existence of subproducts one can simply require that the 
inequality in 4.1 1 hold for all X > 0, Y» A > . 



Geometric criterion 



4.12 



Consider the set of all possible linear combinations 
Vj/lj + v 2 X 2 + V3A3 with X g > . We call it proper cone if it 

lies strictly within one half of the linear space of v ; : 




4.13 

In such a case the integral 4.9 and the product 4.5 are well- 
defined. 

The generalization to more than three factors is straightfor- 
ward. 

Note that the property of a set of vectors to span a proper 
cone is invariant under arbitrary linear transformations. 



Non-linear deformations of coordinates 



4.14 



As a next step towards the general case, notice that exis- 
tence of a distribution is a fact that survives any non-linear 
smooth mapping of coordinates. To rewrite the above criterion 
in a form that is "covariant" under such transformations, con- 
sider a neighbourhood O of the origin and perform a smooth 
( C°° ) deformation of the coordinates so that 



v g p^v g p + 0(p 2 ) = D g (p). 
Fig. 4.13 is modified as follows: 



4.15 




4.16 
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Then Eq.4.5 becomes: 
J dp <p(p) 



supp^C O 



(D { (p) - iO)(D 2 (p) - iO)(D 3 (p) - /0) 



. 4.17 



singular in the sense of differential geometry) within (9, which 
means that 



The singularities of the three factors are now localized on non- 
linear manifolds described by non-linear equations: 

ji g ={p\D g (p) = 0}. 4.18 

The integration measure is transformed with a Jacobian that is 
smooth within O and it is absorbed into the redefined test 
function. The region O is deformed accordingly. 

Moreover, the choice of the origin of coordinates cannot af- 
fect existence of the product, so the origin can be shifted arbi- 
trarily. Denote the point of intersection of the three manifolds 
4.18 as p Q . Since v g = VD g (p Q ) , the criterion 4.11 is rewrit- 
ten as follows: 

A 1 VD 1 (p )+A 2 VD 2 (p )+A 3 VD 3 (p )?£0 foralU g >0. 4.19 

Criterion for any number of factors 4.20 

Consider now a product containing an arbitrary number of 
Sohotsky-type factors: 



1 



jJ Y (D g <p) -/op 



4.21 



Without loss of generality one can assume that there is only 
one point p (not a higher-dimensional manifold) where the 

singular manifolds n g intersect all at once; otherwise by a 

smooth change of coordinates one can deform the intersection 
n y = r\ g< = Y n g into a linear subspace (within (9), and then 

only the coordinates that are transverse to n y will matter. 

By analogy with the preceding examples, one writes down 
the following criterion of existence of the product 4.21 at p : 

^l g VD g (p )*0 forall X g >0, 2X >0 - 422 

gEy gEy 

If all subproducts are exist, then it is sufficient to consider only 
A g >0. 

The heuristic meaning of 4.22 is the same as in the example 
above: as is discussed in 4.27 the vector VD(p ) describes the 

dangerous direction for the distribution (D(p) — iO)~" at its 
singular point p , i.e. the only direction in which its Fourier 

transform does not decrease rapidly, thus jeopardizing conver- 
gence in that direction of Fourier-transformed integrals with 
test functions. 

Recall that the above reasoning allows one to prove the cri- 
terion only in situations where all the non-linear singularities 
can be "flattened" simultaneously by a smooth transformation 
to the form of 4.13. Nevertheless, it is useful to summarize the 
experience gained in a intuitive geometrical form. 

Wave front of a Sohotsky-type distribution 4.23 

Consider a distribution (D(p) — i0)~" defined within a re- 
gion O. Its singularities are localized at the points where 
D{p) = . Such points constitute a manifold, 
Jt = { p\D(p) = 0} . Assume that ji is well-behaved (i.e. is non- 



VD(p)\ D(p)=Q *0 



4.24 



Construct the wave front of the distribution (D(p) — i0)~" by 
attaching the vector VZ)(p ) (describing the "dangerous" di- 
rection; cf.also Sec. 4.27 below) to each singular point p G n : 

VO(Po) 



Po > 4.25 
Now the criterion 4.22 can be reformulated in geometric terms: 



Geometric form of the criterion 



4.26 



For each point p where some of the factors of the product 

4.21 are singular, select all those and only those factors that 
are singular at that point. If the corresponding dangerous di- 
rections VD g (p Q ) span a proper cone — and if this holds true 

for any such point p — then the product is integrable with 

arbitrary test functions localized within O. 

Why VD(p ) is a dangerous direction 4.27 

The above analysis was based on a deformation of coordi- 
nates that made all singular manifolds flat simultaneously. 
This is not always possible (recall osculating singularities, 
Sec. 2.20). Therefore, it is useful to consider the microlocal 
mechanism that ensures existence of the products in terms of 
Fourier transforms without invoking global deformations. It is 
at this point that the notion of wave front becomes genuinely 
non-trivial. 

Consider again the distribution A(p) = (D(p) — iO) a near 
a point p where D(p Q ) = but VZ)(p )^0. It turns out 
that any direction n other than the one described by V£)(p ) 
is not dangerous at p in the following sense: one can find a 
test function ip{p) that is localized sufficiently closely around 
p and such that the Fourier transform 

A^(x) = j dpe" ¥ V(p)A(p) decreases rapidly (i.e. faster than 

any power) in the direction of n, i.e. A^(An) ->0 faster than 

any power X~ N as X -» +oo . 

The proposition is quite intuitive: behaviour at infinity of 
Fourier transforms is determined by the behaviour of the origi- 
nal in an arbitrarily small neighbourhood of the origin. The de- 
nominators of the distributions we deal with are linear func- 
tions up to 0(p 2 ) terms which are — as the preceding exam- 
ples suggest — negligible if one focuses on a sufficiently small 
neighbourhood, so that the distributions resemble the standard 
Sohotsky distribution that depends on the coordinate directed 
along VZ)(p ) . One may also note that multiplication by a 

smooth function does not change this property. 

The formal proof which is a straightforward technical im- 
plementation of the above heuristic idea is given in Sec. 4.28 
below. 

A somewhat stronger version of the above property will be 
useful: the rapid decrease A^ (Aw) -> in A is uniform with re- 
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spect to all directions n outside some (arbitrarily narrow) cone 
C containing WD(p Q ) . This is established as follows. First, 

the proof of the above proposition allows one to obtain bounds 
that are, in fact, uniform with respect to all directions from a 
sufficiently narrow cone C n around n. Then one notices that 

directions can be represented as points of the unit sphere which 
is compact; it follows that for any cone C one can select a fi- 
nite number of directions n t whose cones C n contain all di- 
rections outside C . Choosing the largest bound among those 
corresponding to C n , one obtains the required uniform one. 

Proof of the proposition of Sec. 4.27 4.28 
Let the test functions ip(p) be localized within a small 
neighbourhood O of the point p Q . By a shift of coordinates 
one ensures that p Q = . Consider the integral 



jdpe 



-ilnp 



V(P)- 



1 



4.29 



(L( P )+ior 

Split p = (p ± ,P\\) where the first component is orthogonal to 
the manifold L(p) = while the second, tangential to it at 0. 
The Taylor theorem guarantees that L(p) = p- VL(0) + 0(p 2 ) 
for p -» . Therefore, choosing O small enough allows one to 
find a smooth interpolating function 'P(p) such that 
< P(p) = p - VL(0) for all p sufficiently large, and such that 
W(p) = L(p) within O. Without loss of generality, the fol- 
lowing reasoning assumes L(p) = W(p) . 
There are two cases to be considered: 

A. n is directed at an angle to VL(0) ; 

B. The directions of n and VL(0) are exactly opposite. 
Case A. Let us exhibit the fact that the oscillations of 

the exponential along the singular manifold L(p) = for 
X -> +oo suppress the integral as in the case of a smooth func- 
tion — irrespective of the singular behaviour in the transverse 
directions. To this end, split p ± = (p',p") where p is directed 

along the projection of n onto the subspace of p ± , so that p" 
is orthogonal to both n and VL(0) . In the coordinates 
(p',p",P\\) , one has n = (n',0,«n) . Finally, choose new coordi- 
nates (k',k",k») so that: 



K H = x ¥(p), h'k' = n'p' — W||/c || , k" = p" . 



4.30 



One can see that the closer n to VL(0) , the smaller n' , and 
the closer ip(p) should be to p • VL(0) , which requires a 
smaller O. 

In the new coordinates, the integral becomes 

-il(n..K..+nK') , . 1 



4.31 



The Jacobian of the coordinate transformation was absorbed 
into the test function ip(ic) . The integrand is smooth in k' so 
that one can perform the standard integration-by-parts trick 
(Sec. 3.8) to establish the rapid decrease of the entire integral. 

Case B . Here one makes use of the fact that one deals 
with a "slightly deformed" Sohotsky distribution whose Fourier 
transform is a "slightly deformed" 0-function with essentially 
similar nullification properties. 



Similarly to case A, change the coordinates in 4.29 as fol- 
lows: /C|| = L(p), k ± = p ± . Denote N = —n^ > . Then using 

4.8 (with the integration variable changed to A), rewrite the 
integral as: 

I dAj"ck||Ck ± i/</c)exp!|A/C|| +AM2(/c)-An L /c_ L J , 4.32 

XN 

where Q(/c) = p\(ic) — K n . If the neighbourhood O was chosen 

small enough, one can ensure that |Q'(/c)| < <5 < 1 . This allows 
one to integrate by parts: 

i\AK H +XNQ(K)] 



= J 



1 



A + A/VQ'(k) 



{...}: 



etc. 



4.33 



The new factor does not cause problems because it is smooth 
and the denominator does not turn to zero. In this way one ob- 
tains bounds for the integral over k which, after integration 
over A, make manifest the rapid decrease of the integral. This 
completes the proof. 



Osculating singularities 



4.34 



This is the situation when the singular manifolds 
D g (p) = touch rather than intersect in a general fashion. We 

are going to show that the criterion obtained above works here 
too. 

It is sufficient to consider the simplest case of 2- 
dimensional Euclidean p = (p { ,p 2 ) and the product of two 
factors: 



A 1 (p)A 2 (p), A g (p) = (D g (p)-iO)-\ 
Di(p) = Pi-iO, D 2 (p) = p x - p\ -iO. 
The dangerous directions are shown in 4.36: 



4.35 



Pl = P2 



Pl = 



Pi 



Pl = p 2 =0 



4.36 



It is natural to suppose that since the quadratic terms played no 
role in Sees. 4.14^.20 they should not be important here, too. 

Because the singular manifolds do not intersect away from 
p = , it is sufficient to consider an arbitrarily small neigh- 
bourhood of that point. 

The dangerous direction at p = , the same for both fac- 
tors, is 

VD 1 (0) = VD 2 (0) = (1,0) = V . 4.37 

Take a cone C sin „ of directions around V. Choose it narrow 

enough so that C sing + C sing — which is the cone spanned by 

all possible pairs of vectors from C sing — is a proper cone. 

Furthermore, according to the last paragraph of Sec. 4.27, 
given C sing , one can find a pair of test functions ip\(p) and 

ip 2 ( P) localized around k = such that Fourier transforms 

A (x) of the distributions ip (p)XA (p) are rapidly and 
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uniformly decreasing in all directions except those in the cone Applications 



sing • 



Because the functions tp g (p) are smooth and non-zero in a 

small neighbourhood of p = , it is sufficient to prove exis- 
tence of the product 



X 



I\{p)-iO D 2 (p)-i0 

Consider the integral 

r , , Vi(p) v V>2(p) v , s 
J D { (p)-iO D 2 (p)-i0 



4.38 



4.39 



for an arbitrary test function <p( p) (whose localization proper- 
ties are not important because the distributions themselves are 
localized as needed). Its representation in terms of Fourier 
transforms is as follows: 



Jd 2 xd 2 vA 1 (x)A 2 (y)^(x + y). 



4.40 



One sees that the only directions in the 4-dimensional space of 
the aggregate integration variable (x,y) where the product 

A 1 (jc)A 2 (y) does not exhibit a rapid decrease are the ones 
that lie within the direct product C sing X C sing , i.e. the direc- 
tions of the form (x,y) = AfV^.V^) with V t £ C sing . However, 
the construction of C sing ensures that the integrand in such di- 
rections is suppressed by the rapid decrease of the test func- 
tion. This proves absolute convergence of the integral 4.40 and, 
therefore, existence of the products 4.38 and 4.35. 

The proof immediately generalizes to the case of any num- 
ber of factors. 

To summarize, the criterion of Sec. 4.20 is valid irrespective 
of the pattern of intersection of the singular manifolds Ji g — 

whether they intersect in a transverse fashion or osculate — as 
long as the other conditions are satisfied (most notably, 
smoothness of all manifolds Ji g , or VD g (p) ^ within G). 



Extension to asymptotic expansions 



4.41 



One can immediately extend the above results to obtain a 
criterion of existence of asymptotic expansions (cf. Sec. 3.14). 
It will be sufficient to consider the case of two factors. 

Let k be a parameter that goes to zero, k -* +0 . Consider a 

product of Sohotsky-type distributions (D g (p,ic) - iO) " g . As- 
sume that for all (p,ic) £ Ox (-s,k ) (with some e>0), 
D g (p,K) are smooth in (p,/c) and VD g (p,ic) ^ . Then if 

the product of (D g (p,0) - i0)~ N exists for arbitrarily large N 

in the "natural" sense (i.e. in virtue of the Hormander crite- 
rion) 1 , then the asymptotic expansion in k -» in the sense of 
distributions of the original product is equal to the product of 
the asymptotic expansions of factors (recall that the latter ex- 
pansions are simply the Taylor expansions). 

The proof is a straightforward modification of the reasoning 
in Sec.4.28 in the spirit of Sec. 3.14. 



' which means that the product of Taylor expansions of factors is well- 
defined to all orders in k. 



Landau equations 



5.1 



The original Landau equations [16], [17] provided an an- 
swer to this question: at what values of external parameters — 
masses and momenta — does a given multiloop diagram cease 
to be an analytic function? The problem we are considering 
now is phrased differently: given an integrand of a multiloop 
diagram for some fixed values of external momenta and 
masses, where in the space of integration momenta does the 
integrand (the product of propagators) cease to be a well- 
defined distribution? The close relation of the two questions is 
clearly seen from the point of view of the theory of asymptotic 
operation which establishes a connection between the follow- 
ing two effects: 

(i) non-analyticity of the integrated diagram as a function of 
external parameters at some values of the latter (i.e. a presence 
of, say, logarithmic terms in the corresponding asymptotic ex- 
pansion near those values) and 

(ii) non-integrability of the Taylor-expanded integrand (for 
a discussion see [1], [4]). 



Derivation 



5.2 



Consider a product of propagators of the form 2.2. It need 
not be a complete integrand — just a subproduct whose fac- 
tors' s causal singularities pass over a region O C P . Consider 

the points where the causal singularities of all the factors from 
G overlap. Such points are solutions of the system of equations 



gee 



5.3 



However, at some of those points the product exists, and at 
others, it does not — as described by the criterion of Sec. 4.20. 
Recalling that p is the collection of all loop momenta, 
p = (p l ,p 2 ,...p l ) , and that the gradient operator is then repre- 
sented as 2.13, one rewrites 4.22 as follows: 



dl 2 g ( P ) 

M 

M 2 g (p) 



*0 forall A g >0, 2 A s >0 



5.4 



The points p = (p l ,p 2 ,. .. pj) that satisfy 5.3 and 5.4 are the 

points where the product exists in a natural sense (as explained 
in Sec. 3. 13) — such singularities are traditionally called non- 
pinched , while p satisfying 5.3 but not 5.4 constitute a mani- 
fold in the integration space where the product does not exist 
— the manifold of pinched singularities. 

The system of Landau equations is obtained if one replaces 
inequalities by equalities in 5.4, and reinterpretes it as a de- 
scription of pinched singularities rather than conditions for ex- 
istence of the product. 



Non-covariant gauges 



5.5 



In the theory of non-covariant light-cone gauges (nA = 

with n 2 = ) one encounters propagators of the following 
form: 
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p 2 +i0 



nuPv +n vPu 



3 fiV 



n p 



5.6 



where n 2 = 0, n° > 



An important technical problem (cf. 

e.g. [31]) is the singular structure of Feynman diagrams in- 
volving such propagators in view of the additional singularity 
in 5.6. The so-called Mandelstam-Leibbrandt prescription (for 
a review see [34]) consists in using 

,0 



np^np + iOX sign p 



5.7 



Using the formalism developed above, it is not difficult to see 
why this prescription works: 

Indeed, the singular manifolds corresponding to the two 
factors in the denominator of 5.6 are described by p 2 = and 
np = 0, and intersect over the line p^n . Consider any point 
p on that line and evaluate the gradients: 



V (signp Xrcp) = signp Xn 



5.8 



(The overall sign is switched so as to make the imaginary in- 
finitesimal parts the same sign, as required by the formalism of 
the preceding sections.) 

One sees that the two gradients are the same up to a posi- 
tive coefficient — irrespective of whether p is on the upper or 
lower half-cone. The criterion 4.26 then tells us that the prod- 
uct 5.6 with the prescription 5.7 is a well-defined distribution. 
This also means that the wave front of the propagator 5.6 on 
the light cone is not different from that of the simple scalar 
propagator. From the point of view of studying intersections of 
causal singularities of several such propagators, this means that 
there will be no new pinches etc. as compared with the covari- 
ant case. Of course, the propagator 5.6 has an additional sin- 
gularity np = outside the light cone. But since such addi- 
tional singularities are localized on a linear manifold, studying 
their intersections presents no difficulty. 

A reasoning of this sort allows one to largely avoid explicit 
calculations when studying the singular structure of diagrams 
(cf. examples of such calculations in [31]). 



Quasi-Euclidean asymptotic regimes 



5.9 



As a last application, we consider the class of asymptotic 
regimes which, although not allowing Wick rotation, never- 



theless result in the same formulas as in the Euclidean case. A 
simple example of such a regime was given in Sec. 2.21. 

It is not difficult to understand (cf. Fig. 2.24) that no modi- 
fications (no new counterterms) to the formulas of Euclidean 
asymptotic operation are necessary as long as the following 
holds true in the asymptotic limit corresponding to a particular 
regime: Apices of light cones corresponding to various propa- 
gators should either merge in the asymptotic limit, or stay away 
from each other at a non-light-like separation; they are also not 
allowed to approach mass shell surfaces in the asymptotic 
limit. More precisely, if, as in [4], m, q and M, Q are small 
and large masses and external momenta of the problem, re- 
spectively (for definiteness, all the momenta enter the dia- 
grams), and if one formulates the asymptotic regime in terms 
of m,g-»fj with M,Q fixed, then q should vanish compo- 
nentwise whereas Q are allowed to be arbitrary non-lightlike 
provided the asymptotic configuration does not correspond to 
any of the thresholds. More precisely, in the context of a con- 
crete 1PI diagram G, let C a be a cut of the diagram with Q a 

being the momentum entering any one of the two resulting 
components of G, and let M a be the sum of masses of the 

lines cut by C a ; further suppose that Q 2 ^ M 2 for any cut a 
in the asymptotic limit ( m = q = ). Then the asymptotic ex- 
pansion is given by the same formulas as in the case of purely 
Euclidean asymptotic regimes [4], [5]. 

Finally, it should be emphasized that, strictly speaking, the 
pinch/non-pinched classification of the singularities can be in- 
terpreted in the context of asymptotic operation as equality to 
zero of counterterms at the corresponding points of singular 
manifolds. Therefore, one might proceed with the construction 
of non-Euclidean asymptotic operation ignoring such issues, 
and then rediscover nullification of counterterms for non- 
pinched singularities by an inspection of the resulting integrals. 
As a matter of fact, analytical study of the most complex types 
of non-Euclidean singularities (boundary points of pinched 
manifolds) may be psychologically easier if one ignores the 
pinch problem altogether. But the existing treatments start 
with the pinch/non-pinched classification (cf. [35]), so it would 
have been inconvenient to postpone establishing a connection 
with what is already known until the new theory is completely 
developed. 
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